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ABSTRACT 

This  is  the  second  of  two  papers  examining  Spitzer  Infrared  Spectrograph  observations  of  the  ultraluminous 
X-ray  source  (ULX)  in  Holmberg  IT  Here,  we  perform  detailed  photoionization  modeling  of  the  infrared 
lines.  Our  analysis  suggests  that  the  luminosity  and  morphology  of  the  [Oiv]  25.89  /xm  emission  line  is 
consistent  with  photoionization  by  the  soft  X-ray  and  far  ultraviolet  (FUV)  radiation  from  the  accretion  disk 
of  the  binary  system  and  inconsistent  with  narrow  beaming.  We  show  that  the  emission  nebula  is  matter- 
bounded  both  in  the  line-of-sight  direction  and  to  the  east,  and  probably  radiation-bounded  to  the  west.  A 
bolometric  luminosity  in  excess  of  10"^^  erg  s“^  would  be  needed  to  produce  the  measured  [Oiv]  flux.  We 
use  modeling  and  previously  published  studies  to  conclude  that  shocks  likely  contribute  very  little,  if  at  all, 
to  the  high-ionization  line  fluxes  observed  in  the  Holmberg  II  ULX.  Additionally,  we  And  that  the  spectral 
type  of  the  companion  star  has  a  surprisingly  strong  effect  on  the  predicted  strength  of  the  [Oiv]  emission. 
This  flnding  could  explain  the  origin  of  [Oiv]  in  some  starburst  systems  containing  black  hole  binaries. 
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1.  INTRODUCTION 

Despite  over  two  decades  of  research  (Fabbiano  1989),  lit¬ 
tle  is  known  about  the  intrinsic  properties  of  ultraluminous 
X-ray  sources  (ULXs).  They  are  characterized  by  extremely 
high  X-ray  luminosities,  some  reaching  10"^^  erg  s“^  These  lu¬ 
minosities  cannot  be  easily  explained  by  normal  stellar-mass 
black  hole  (BH)  systems,  such  as  those  found  in  our  Galaxy. 
Beamed  or  super-Eddington  stellar-mass  systems  have  been 
proposed  as  possible  models.  Alternatively,  ULXs  can  harbor 
intermediate-mass  black  holes  (IMBHs).  For  a  complete  discus¬ 
sion  of  various  ULX  mechanisms,  see  for  example  the  review  pa¬ 
pers  of  Miller  &  Colbert  (2004),  Soria  et  al.  (2005),  and  Roberts 
(2007).  Some  ULX  X-ray  spectra  show  distinct  features,  which 
have  been  sometimes  interpreted  as  evidence  for  the  IMBH  sce¬ 
nario  (Miller  et  al.  2004).  However,  the  super-Eddington  mecha¬ 
nism  (maybe  together  with  mild  beaming)  has  gained  ground  in 
more  recent  years,  backed  both  by  theoretical  and  observational 
evidence  (Freeland  et  al.  2006;  Done  &  Kubota  2006;  Kuncic 
et  al.  2007;  King  2008;  Soria  &  Kuncic  2008;  Berghea  et  al. 
2008). 

While  ULXs  are  easy  to  detect  in  the  X-rays,  they  are  often 
faint  or  undetected  in  other  energy  bands,  making  them  very  dif- 
flcult  to  characterize.  The  ionization  nebulae  surrounding  some 
ULXs  have  become  critical  for  understanding  the  properties 
of  the  central  source,  such  as  the  intrinsic  luminosity.  For  in¬ 
stance,  Kaaret  et  al.  (2004,  KWZ  hereafter)  used  the  optical 
Hen  recombination  line  (A.4686)  to  characterize  the  intrinsic 
X-ray  luminosities  of  ULXs,  as  was  previously  done  for  a  “nor¬ 
mal”  X-ray  binary  by  Pakull  &  Angebault  (1986).  Even  [Ne  v] 
A.3426,  a  much  higher  excitation  line  has  been  recently  detected 
in  a  ULX  (Kaaret  &  Corbel  2009).  However,  we  note  that  these 
and  analogous  studies  are  limited  to  the  optical  regime. 

The  very  luminous  ULX  in  the  dwarf  galaxy  Holmberg  II 
is  one  such  object  that  has  been  studied  extensively  in  both 
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the  X-ray  and  optical  bands  (Zezas  et  al.  1999;  Lehmann  et  al. 
2005;  Goad  et  al.  2006;  Abolmasov  et  al.  2007).  It  is  located 
inside  an  ionized  nebula  (the  “Eoot  nebula”),  and  shows  high- 
ionization  optical  emission  lines  coincident  with  the  X-ray 
source  (Pakull  &  Mirioni  2002;  Abolmasov  et  al.  2007).  The 
nebula  appears  to  be  perturbed  by  a  moderate  velocity  held  of 
~50  km  s“^  (Lehmann  et  al.  2005).  The  optical  counterpart  of 
the  ULX  is  a  stellar-like  bright  source  located  in  the  “Heel” 
of  the  nebula,  having  emission  consistent  with  a  B  supergiant  to 
a  main- sequence  O  star  (Kaaret  et  al.  2004).  Morever,  extended 
radio  emission  was  detected,  approximately  coincident  with 
the  ionized  nebula  (Miller  et  al.  2005).  The  high-ionization 
recombination  line  Hen  A.4686  suggests  photoionization  by 
the  ULX  (KWZ),  and  the  radio  emission  is  inconsistent  with 
emission  from  a  supernova  remenant  (SNR)  or  an  H  ii  region. 

In  our  first  paper  (Berghea  et  al.  2009,  Paper  I  hereafter)  we 
presented  Spitzer  Infrared  Spectrograph  (IRS)  observations  of 
the  Eoot  nebula  in  Holmberg  II .  The  low-resolution  spectral 
maps  presented  in  that  paper  suggest  that  the  high-ionization 
[O  iv]  line  emission  is  coincident  with  the  ULX  position.  Indeed 
we  found  the  mid-IR  line  ratios  for  this  object  to  be  very  similar 
to  those  seen  in  standard  active  galactic  nucleus  (AGN;  Genzel 
et  al.  1998;  Sturm  et  al.  2002;  Dale  et  al.  2006).  Our  preliminary 
analysis  suggested  that  [O  iv],  like  the  He  ii  recombination  line, 
can  be  used  to  constrain  the  spectral  energy  distribution  (SED) 
of  the  ULX  in  both  the  UV  and  X-rays,  thereby  constraining  the 
intrinsic  X-ray  luminosity.  To  this  end,  we  presented  archival 
X-ray,  ultraviolet  (UV),  and  infrared  (IR)  data  and  constructed 
a  detailed  SED  of  the  ULX  from  which  photoionization  models 
can  be  used  to  break  the  degeneracy  in  X-ray  models  and  extract 
a  reliable  bolometric  luminosity  for  the  source. 

In  this  second  paper,  we  perform  detailed  photoionization 
modeling  for  the  Holmberg  II  ULX.  We  use  the  SED  obtained 
in  Paper  I  as  the  ionizing  source,  and  compare  the  predictions 
with  the  observed  IR  Spitzer  spectrum.  Section  2  presents  a 
summary  of  the  earlier  results,  which  are  then  used  in  Section  3 
to  perform  the  photoionization  modeling.  To  better  constrain  the 
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Table  1 

X-ray  Model  Fits 


Model 

(1) 

N// 

(2) 

(10^0  cm-3) 

k7in/r  1 

(3) 

(keV) 

r/r2 

(4) 

Ax^/dof 

(5) 

L 

(6) 

(10^*^  erg  s“^) 

MCD  flux 

(7) 

Modified  PLMCD 

3.26  ±  0.56 

0.38  ±  0.02 

2.54  ±  0.02 

1.02/135,0.974/713 

1.11 

0.49 

PLMCD 

10.4  ±  0.70 

0.27  ±  0.02 

2.42  ±  0.05 

1.1386/849 

2.28 

0.12 

BPL 

5.49  ±  0.78 

1.73  ±0.10 

2.57  ±  0.02 

1.0316/849 

1.41 

0.60 

Notes.  (1)  X-ray  model,  as  described  in  Section  2.  (2)  Intrinsic  hydrogen  column  density.  The  Galactic  column  density  from 
KWZ  (3.42  X  10^^  cm“^)  was  added  separately.  (3)  Inner  disk  temperature  for  the  MCD  component.  For  the  BPL,  this  is 
the  first  (low  energies)  photon  index  parameter.  The  break  is  at  1.0  keV.  (4)  Photon  index  for  the  power-law  component.  For 
the  broken  power-law  model  this  is  the  second  (high  energies)  photon  index  parameter.  (5)  Reduced  values  for  the  fit  and  the 
number  of  degrees  of  freedom.  For  the  modified  PLMCD  model,  we  quote  two  values  for  each  separately  fitted  component  (see  the  text 
for  details).  (6)  Unabsorbed  (intrinsic)  luminosities  between  0.1  and  10  keV.  (7)  MCD  component  unabsorbed  flux  as  fraction  of  the 
total  flux  for  the  PLMCD  model.  For  the  other  two  models,  this  is  simply  the  ratio  of  the  flux  between  0.1  and  1.0  keV  to  the  flux  in  the 
whole  range  (0.1-10  keV). 


physical  properties,  we  investigate  how  the  predictions  change 
when  we  vary  the  parameters:  the  accretion  disk  temperature, 
gas  density,  metallicity,  and  the  stellar  companion.  We  try  both 
the  radiation-  and  matter-bounded  geometries  to  explain  the 
high-  and  low-ionization  IR  lines.  Finally,  in  Section  4,  we 
explore  shock  models  as  a  possible  alternative  to  ionization.  In 
Sections  5  and  6,  we  discuss  our  results  and  we  present  our 
conclusion. 

2.  SUMMARY  OF  OBSERVATIONAF  RESUFTS 

The  study  of  UFXs  is  limited  almost  exclusively  to  the  X-ray 
band,  where  they  appear  very  bright  and  uncontaminated.  Stud¬ 
ies  of  these  objects  at  other  wavelengths  are  difficult  because 
the  emission  is  generally  either  undetected  or  contaminated  by 
circumstellar  or  circumnebular  emission,  making  an  estimate 
for  the  intrinsic  bolometric  luminosity  very  difficult.  This  de¬ 
ficiency,  especially  in  the  UV,  results  in  degeneracy  in  X-ray 
models  of  the  accretion  process  that  cannot  be  broken.  Our 
analysis  of  the  X-ray  emission  from  the  UFX  in  Holmberg  II 
is  an  example  of  such  degeneracy.  In  Paper  I,  three  physically 
distinct  models  were  presented  to  fit  the  X-ray  spectra  for  the 
Holmberg  II  UFX.  The  resultant  parameters  from  the  fits  are 
given  in  Table  1.  Figures  1-3  show  the  spectra  and  the  fitted 
components  for  all  three  models.  A  brief  summary  of  the  three 
models  is  as  follows. 

1.  Power-law  plus  multicolor  disk  (PLMCD)  model.  This  is  a 
two-component  model  widely  used  to  describe  the  emission 
from  accreting  BHs.  As  the  name  suggests,  the  model 
consists  of  a  multicolor  disk  component  (MCD)  and  a 
power-law  component  for  the  hard  tail.  As  discussed  in 
Paper  I,  when  extrapolated  to  softer  energies  (the  far  UV) 
this  model  falls  apart  since  the  power  law  dominates  in  the 
UV.  Indeed  the  model  is  unrealistic  since  it  is  physically 
impossible  to  have  the  X-ray  photons  Comptonized  to 
lower  energies.  This  model  provides  a  good  fit  to  the 
X-ray  spectrum  and  is  used  extensively  in  other  X-ray 
studies.  However  we  emphasize  that  it  is  not  physically 
meaningful  (e.g..  Done  &  Kubota  2006).  The  PFMCD 
model  is  shown  here  as  a  base  from  which  modified,  but 
more  physically  meaningful  models  are  derived. 

2.  The  modified  PLMCD  model.  The  modified  PFMCD  model 
is  also  a  two-component  model,  but  one  for  which  the  far- 
UV  extrapolation  is  from  the  disk  component  rather  than  the 
power-law  component.  This  yields  a  much  more  physically 
meaningful  interpretation  for  the  emission. 


PLMCD  Model 


Figure  1.  PLMCD  model  fit  to  the  XMM-Newton  spectrum.  This  is  the  typical 
model  used  for  ULXs:  a  MCD  plus  a  power-law  at  high  energies. 


3.  Broken  power-law  (BPL)  model.  The  spectra  were  also 
fit  with  a  broken  power-law  model,  which  is  much  more 
flexible  when  extrapolated  to  the  UV.  Such  a  model  would 
be  produced  if  jet  emission  dominated  the  full  SED.  We  note 
that  there  is  currently  no  direct  evidence  for  jets  in  UFXs.  It 
is  interesting  however  that  one  of  the  few  sources  associated 
with  radio  emission,  the  famous  UFX  in  NGC  5408,  was 
shown  by  Kaaret  et  al.  (2003)  to  be  well  fitted  by  a  broken 
PF  (see  also  Soria  et  al.  2006;  Kaaret  &  Corbel  2009).  For 
more  details  on  our  motivation  for  these  models  we  refer 
our  reader  to  Paper  I. 

As  can  be  seen  from  Table  1,  all  three  models  provide  rea¬ 
sonably  good  fits  to  the  X-ray  data.  This  degeneracy  can  only 
be  broken  through  comprehensive  knowledge  of  the  shape  of 
the  UV  continuum  and  detailed  photoionization  modeling  of 
high-excitation  lines.  In  the  former  case,  the  three  X-ray  mod¬ 
els  can  be  extrapolated  to  the  UV  and  then  compared  with 
observational  data.  The  best-fit  extrapolated  model  to  the  UV 
spectrum  provides  some  indication  of  the  appropriate  physical 
accretion  scenario  for  the  observational  data.  However,  con- 
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Modified  PLMCD  Model 


Figure  2.  Modified  PLMCD  model  fit.  Here  the  MCD  component  and  the  power 
law  are  fitted  separately  below  and  above  1  keV,  respectively. 


tamination  to  the  UV  emission  by  circumstellar  or  circumneb- 
ular  continuum  precludes  a  pure  fit.  Thus,  in  addition  to  this 
photometric  data,  high-ionization  spectral  lines,  which  require 
high-energy  UV  photons  to  ionize,  can  be  used  in  combination 
with  photoionization  modeling  to  check  the  robustness  of  the 
fit.  In  Paper  I,  we  also  presented  multi- wavelength  photometric 
data  from  archival  Spitzer,  XMM-Newton,  and  GALEX  observa¬ 
tions  of  the  Holmberg  II ULX,  and  also  Hubble  Space  Telescope 
(HST)  and  Very  Large  Array  photometric  data  from  the  literature. 
These  data  were  used  to  construct  a  detailed  SED  of  the  ULX 
(Eigure  4).  All  three  X-ray  models  enumerated  above  are  shown 
in  the  plot  as  dotted,  thick  solid,  and  dashed  black  lines,  re¬ 
spectively.  The  thick  solid  black  line  in  Eigure  4  represents  the 
“base  SED  model”  which  is  the  input  spectrum  for  the  prelimi¬ 
nary  photoionization  modeling  explored  below.  As  discussed  in 
Paper  I  and  shown  in  Eigure  4,  the  modified  PLMCD  model  is 
chosen  as  the  X-ray  model  for  the  base  model,  since  it  is  the 
only  extrapolated  model  that  fits  the  optical  and  UV  data  well. 
Tables  1  and  2  provide  the  X-ray  spectral  fitting  results  and  the 
IR  line  measurements  from  the  first  paper  that  will  be  the  basis 
of  the  modeling  in  the  current  publication. 

3.  CLOUDY  MODELING  OE  THE  IR  LINES 

To  perform  detailed  photoionization  modeling  with 
CLOUDY,^  a  number  of  input  parameters  need  to  be  defined. 
Specifically,  the  geometry  of  the  cloud  needs  to  be  described, 
the  density  and  metallicity  of  the  gas  needs  to  be  specified,  and 
an  input  SED  needs  to  be  established.  The  following  outlines 
the  specific  input  parameters  used  for  the  preliminary  CLOUDY 
run.  These  parameters  are  later  varied  to  check  the  final  results. 

1.  Gas  cloud  geometry.  The  input  geometry  must  describe  the 
shape  and  thickness  of  the  gas  cloud.  Shape  includes,  for  ex¬ 
ample,  a  spherical,  cylindrical,  or  plane-parallel  geometry. 
The  cloud  thickness  defines  a  radiation-  or  matter-bounded 


^  Photoionization  modeling  was  performed  using  ver.  07.02.01  of  CLOUDY, 
described  by  Ferland  et  al.  (1998). 


BPL  Model 


Figure  3.  Broken  power-law  model  fit.  The  break  is  at  1  keV,  the  same  energy 
where  the  two  components  are  separated  in  the  modified  PLMCD  model  in 
Figure  2. 

geometry.  In  their  photoionization  modeling,  KWZ  as¬ 
sumed  a  spherical  geometry  based  on  resolved  optical 
lines,  and  assumed  the  nebula  is  radiation-bounded.  How¬ 
ever,  they  noted  that  the  lower-ionization  line  [O  i]  is  only 
detected  west  of  the  ULX,  thereby  suggesting  a  matter- 
bounded  nebula  in  other  directions.  Though  we  also  de¬ 
tected  low-ionization  lines  in  the  IR  spectrum  (Table  2  in 
this  paper  and  Eigures  1  and  2  in  Paper  I),  we  begin  with 
the  spherical,  radiation-bounded  geometry  in  an  attempt  to 
reproduce  the  high-excitation  [O  iv]  emission.  We  assume 
the  inner  radius  of  the  cloud  to  be  0.1  pc,  thus  maximizing 
the  ionizing  fiux.  To  constrain  the  size  of  the  ionized  neb¬ 
ula  for  the  radiation-bounded  geometry,  CLOUDY  was  run 
down  to  a  temperature  of  4000  K.  At  this  temperature,  the 
cloud  thickness  is  140  pc,  and  the  column  density  is  4.2  x 
10^^  cm“^,  more  than  10  times  the  X-ray  estimated  value. 

2.  Density.  Density  can  be  specified  as  varying,  uniform  or 
clumpy.  In  the  preliminary  CLOUDY  model,  we  assume  a 
constant  gas  density  of  10  cm“^  and  a  filling  factor  of  1  as 
per  KWZ.  This  density  (and  filling  factor)  was  estimated 
by  KWZ  using  the  optical  line  surface  brightness  profiles 
predicted  by  CLOUDY. 

3.  Metallicity.  Metallicity  is  expressed  in  terms  of  solar 
metallicities.  We  chose  an  input  metallicity  of  0. 1  Zq,  close 
to  that  of  KWZ  and  Pilyugin  et  al.  (2004). 

4.  SED.  The  input  SED  must  specify  the  companion  star  and 
the  X-ray  model  used  to  fit  the  UV  to  X-ray  continuum.  The 
input  SED  used  in  our  preliminary  analysis  is  described  in 
Section  2  and  results  from  the  multi-wavelength  photomet¬ 
ric  and  spectroscopic  observations  presented  in  Paper  I. 
This  so-called  “base  SED  model”  utilizes  the  modified 
PLMCD  model,  extrapolated  to  the  UV.  Eor  the  input  SED 
we  used  a  B2Ib  supergiant  spectrum  to  represent  the  stel¬ 
lar  companion.  The  B2Ib  supergiant  was  chosen  because  it 
better  fits  our  UV  data  than  the  05V  (both  are  limits  to  the 
stellar  spectral  range  suggested  by  KWZ).  The  bolometric 
luminosity  of  this  model  is  1.34  x  10"^^  erg  s“^. 
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Figure  4.  Input  SED  (base  model)  constructed  in  Paper  I  and  used  as  ionizing  source  in  CLOUDY.  The  radio  data  are  from  Miller  et  al.  (2005),  the  U-band  magnitude 
for  the  optical  counterpart  of  the  ULX  is  quoted  from  KWZ.  Our  measured  fluxes  in  Paper  I  using  IRAC,  MIPS,  OM,  and  GALEX  are  upper  bounds,  and  therefore  our 
model  is  below  these  data  points.  The  two  stellar  spectra  correspond  to  the  spectral  types  range  consistent  with  the  colors  and  magnitudes  in  KWZ.  The  three  models 
fit  to  the  X-ray  data  are  shown,  extrapolated  to  the  [O  iv]  edge  at  54.93  eV  (shown  as  a  vertical  line).  The  BPL  model  was  used  to  obtain  a  family  of  broken  power-law 
models,  all  compatible  with  the  X-ray  data  within  the  errors.  We  show  the  lower  and  upper  limits  of  this  family  as  the  lower  and  upper  dashed  lines,  corresponding  to 
r  =  1.3  and  E  =  2.2  in  the  low  energy  slope,  respectively.  The  middle  dashed  line  represents  the  best-fit  BPL  model  shown  in  Table  1.  The  SED  is  constructed  using 
the  preferred  X-ray  model,  modified  PLMCD. 


Table  2 

Measured  Infrared  Lines 


Line 

(1) 

Standard  Aperture 

Small  Aperture 

Flux 

(2) 

S/N  Ratio 
(3) 

L 

(4) 

Flux 

(5) 

S/N  Ratio 
(6) 

L 

(7) 

[Neii]  12.81  fim 

<6.08 

<0.68 

<1.98 

<0.22 

[Neiii]  15.56  /xm 

12.15  ±  3.4 

4.60 

1.35  ±  0.39 

5.46  ±  1.37 

7.6 

0.61  ±  0.15 

[Sill]  18.71  /xm 

2.64  ±  1.14 

4.64 

0.29  ±  0.13 

3.02  ±  1.16 

4.7 

0.34  ±  0.13 

[Nev]  24.32  /xm 

<3.88 

<0.43 

<1.16 

<0.13 

[Oiv]  25.89  /xm 

7.01  ±  1.60 

10.16 

0.78  ±  0.2 

3.63  ±  0.91 

44.4 

0.40  ±  0.1 

[Sill]  33.48  /xm 

<4.39 

<0.49 

<1.53 

<0.17 

[Sin]  34.82  /xm 

21.91  ±  5.48 

9.90 

2.44  ±  0.61 

5.87  ±  1.47 

7.2 

0.65  ±  0.16 

Notes.  We  show  line  fluxes  for  the  standard  (slightly  smaller  than  the  SH  map)  and  the  small  aperture  (4  LH  pixels)  as  defined  in 
Section  2  of  Paper  I.  Fluxes  are  in  10“^^  W  cm“^,  luminosities  (L)  in  10^^  erg  s“^  If  the  measurement  errors  are  smaller  than 
the  absolute  calibration  accuracy  of  25%,  the  latter  were  used.  For  nondetections  we  show  3a  upper  limits. 


3.1.  Results  from  the  Base  Model 

The  results  of  the  first  CLOUDY  run,  with  the  base  models  as 
the  ionizing  source,  are  given  in  Figure  5  (top-left  panel),  as  the 
ratio  of  the  predictions  to  the  observed  values.  The  CLOUDY 
model  predicted  values  for  the  mid-IR  lines  are  close  to  the 
observed  values  for  the  high-ionization  lines.  A  simple  change 
to  the  base  model  such  as  a  slightly  younger  companion  star 
would  be  enough  to  fit  the  [O  iv]  line  observed,  as  the  analysis 
below  indicates.  We  also  show  predictions  for  the  Hen  fiux 
measured  by  KWZ  (2.7  x  10^^  erg  s“^). 

The  low-ionization  lines  are  the  most  discrepant.  This  dispar¬ 
ity  between  the  model  and  the  observations  is  likely  caused  by  a 
matter-bounded  geometry  in  at  least  some  directions  of  the  neb¬ 


ula,  as  already  suggested  by  KWZ.  However,  many  of  the  model 
parameters  are  not  well  known  and  we  explore  how  changes  to 
these  can  effect  the  line-fiux  predictions  in  the  sections  to  follow. 

We  note  that  sensitivity  and  aperture  effects  could  explain  the 
discrepancy  between  the  low-ionization  line  fiuxes  predicted  by 
the  radiation-bounded  base  SED  model  and  the  observed  values. 
Using  the  CLOUDY  output  we  calculated  surface  brightness 
profiles  expected  for  the  mid-IR  lines.  CLOUDY  predicts  a 
spatial  extent  of  20-25  pc  for  the  He  ii  which  is  consistent  with 
the  findings  of  KWZ  (15-26  pc).  The  [Oiv]  emission  from  the 
ULX  should  be  extended  to  a  radius  of  about  ~20  pc  according 
to  the  simulation.  However,  due  to  the  poor  resolution  of  our 
spectral  maps  (one  LH  pixel  =  67.5  pc),  it  is  unclear  whether 
the  model  accurately  predicts  the  spatial  extent  of  the  [O  iv] . 
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Figure  5.  We  show  predictions  for  the  IR  lines  from  the  CLOUDY  modeling  for  various  ionizing  SEDs.  We  also  included  the  Hen  A. 4686  line  measured  by  KWZ 
(2.7  X  10^^  erg  s“^).  They  are  presented  as  model-to-data  flux  ratios.  For  non-detected  lines,  the  plotted  values  are  lower  limits  and  are  marked  with  an  arrow  pointed 
upward.  The  predictions  for  the  [Ne  iii]  and  [Si  ii]  lines  are  <0.01  in  the  BPL  model  (lower-left)  and  are  not  shown.  The  perfect  match  line  (ratio  =  1)  is  shown  on  all 
plots.  The  errors  are  small,  similar  in  size  to  the  data  points,  therefore  are  not  shown.  The  measured  line  fluxes  are  presented  in  Table  2. 


In  addition,  the  CLOUDY  models  suggest  that  while  sensitivity 
may  have  some  effect  on  the  observed  low-ionization  line  fluxes, 
it  is  not  solely  responsible  for  the  high  predicted-to-observed 
line  ratios  for  the  radiation-bounded  base  SED  model. 

3.2.  Varying  the  X-ray  Source  Spectrum  to  the  BPL  and 
PLMCD  Models 

We  next  tried  changing  the  X-ray  part  of  the  source  spectrum. 
We  first  replace  the  modified  PLMCD  model  with  the  PLMCD 
model  described  in  Section  2,  leaving  the  other  CLOUDY 
parameters  the  same.  The  model  is  extrapolated  to  the  edge 
of  the  [Oiv]  line  where  we  introduce  a  break,  so  that  it 
does  not  dominate  the  optical  fiux  from  the  B2  star.  Results 
from  this  CLOUDY  simulation  are  shown  in  the  upper-right 
panel  of  Eigure  5.  The  values  predicted  by  this  model  are 
clearly  not  consistent  with  the  measured  values.  This  model 
predicts  a  strong  [Nev]  line,  which  is  not  detected  in  our 
Spitzer  observations,  and  a  He  ii  luminosity  27  times  larger  than 
measured  by  KWZ.  The  [O  iv]  fiux  predicted  from  this  model 
is  also  an  order  of  magnitude  larger  than  the  observed  fiux.  We 
therefore  rule  out  this  model  as  a  good  estimate  for  the  X-ray 
source  spectrum.  We  note  that  the  PLMCD  model  has  been  used 
by  many  researchers  to  estimate  disk  temperatures  and  even  BH 
masses.  Our  findings  imply  that  it  is  likely  not  an  appropriate 
model  for  ULX  high-energy  spectra. 

We  then  tried  replacing  the  modified  PLMCD  model  with 
the  BPL  model  described  in  Section  2.  Here,  we  found  results 
similar  to  the  PLMCD  model,  with  [Nev]  predicted  to  be 
detected,  and  all  the  observed  lines  well  overpredicted.  We 
therefore  conclude  that  neither  the  PLMCD  model  nor  the  BPL 
model  with  a  radiation-bounded  geometry  provide  a  good  fit  to 


the  observed  IR  spectra.  This  was  expected  because  both  these 
models  predict  a  very  strong  UV  ionizing  fiux. 

3.3.  Varying  Other  Parameters 

We  investigate  here  in  more  detail  how  the  strong  radiation 
from  the  accretion  disk  photoionizes  the  surrounding  interstellar 
medium  (ISM)  and  produces  the  observed  lines.  Starting  again 
with  the  base  SED  model  in  the  radiation-bounded  geometry, 
we  vary  the  input  parameters  to  see  how  they  affect  the  emitted 
spectrum.  We  first  examine  the  dependence  of  the  simulated 
line  luminosities  on  the  X-ray  spectrum  by  changing  the  MCD 
disk  temperature  (kTin),  then  the  gas  density,  metallicity,  and 
the  inner  radius  of  the  cloud.  Linally,  we  explore  how  the 
SED  of  the  stellar  companion  affects  the  line  luminosities. 
The  models  presented  are  based  on  the  multi-wavelength  data 
of  the  Holmberg  II  ULX  and  the  known  properties  of  the 
surrounding  ISM,  but  this  analysis  can  be  applied  to  ULXs  in 
general. 

Varying  the  disk  temperature.  The  disk  temperature  is  first 
altered  by  shifting  the  X-ray  base  SED  model  to  lower  photon 
energies,  such  that  the  disk  luminosity  remains  constant  across 
the  band.  As  can  be  seen  from  Eigure  6(a),  when  the  disk 
component  peaks  at  lower  energies,  the  mid-IR  emission  lines 
become  stronger.  Lor  cooler  disks  (0.1  keV)  the  [Nev]  line 
should  be  detected  in  our  observations.  In  contrast,  the  lower 
ionization  lines  are  not  dependent  on  the  disk  temperature  over 
the  range  usually  found  in  ULXs  (0. 1-0.5  keV).  We  conclude 
from  this  that  varying  the  disk  temperature  will  not  yield  better 
results. 

Varying  the  gas  density  of  the  cloud.  The  results  obtained  by 
varying  the  density  of  the  cloud  are  plotted  in  Eigure  6(b).  The 
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Figure  6.  Dependence  of  the  simulated  lines  from  CLOUDY  on  (a)  MCD  inner  disk  temperature  and  (b)  gas  density  of  the  cloud.  The  vertical  lines  mark  the  base 
model  values:  kTi^  =  0.38  keV,  density  =  10  cm“^.  The  line  colors  are  the  same  but  only  labeled  in  the  left  panel. 


Figure  7.  Dependence  of  the  simulated  lines  from  CLOUDY  on  (a)  cloud  metallicity  and  (b)  the  inner  radius  of  the  cloud.  The  vertical  line  marks  the  base  model 
value  of  Z  =  O.1Z0.  The  line  colors  are  the  same  but  only  labeled  in  the  left  panel. 


line  luminosities  are  not  particularly  sensitive  to  gas  density 
(see  also  KWZ),  though  the  [O  iv]  luminosity  increases  slightly 
as  the  gas  approaches  the  critical  density.  KWZ  found  that  the 
Hen  line  is  clearly  extended  to  15-25  pc  and  consistent  with  a 
density  of  "^10  cm“^.  We  therefore  conclude  that  varying  the 
density  of  the  gas  will  not  do  much  to  improve  the  model. 

Varying  the  metallicity.  The  dependence  of  the  line  lumi¬ 
nosities  on  metallicity  is  shown  in  Figure  7(a).  The  [Oiv]  and 
[Ne  v]  line  dependence  is  linear.  Interestingly,  the  base  model 
with  a  higher  input  metallicity  reproduces  the  measured  [O  iv] 
luminosity  very  well.  In  fact,  a  metallicity  of  0.19  Zq  (oxy¬ 
gen  abundance  of  9.3  x  10“^  in  absolute  value)  reproduces  the 


[O  iv]  to  within  the  calibration  error  of  the  instrument,  but  with¬ 
out  increasing  the  He  ii,  which  is  already  overpredicted  in  the 
current  version  of  the  model.  We  conclude  from  this  that  the 
input  metallicity  may  be  lower  than  it  should  be  (see  Pilyugin 
et  al.  2004).  However,  we  also  note  that  a  higher  metallicity  will 
predict  higher  fluxes  for  the  low-ionization  lines,  resulting  in 
an  even  greater  discrepancy  for  these  lines.  We  also  explored 
the  effects  of  gas  phase  depletion  on  metallicity  and  found  that 
depletion  had  very  little  effect  on  either  the  metallicity  or  the 
resultant  luminosity  of  the  [O  iv]. 

Varying  the  inner  radius  of  the  ionizing  cloud.  Figure  7(b) 
shows  that  the  ionizing  power  of  the  incident  flux  is  diluted  if 
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Figure  8.  Comparison  between  the  base  model,  with  the  B2Ib  star  (continuous  line)  and  the  model  with  the  05 V  star  (dashed  line),  (a)  The  oxygen  ionic  structure  in 
the  cloud  (O^"^  and  O^"^).  (b)  Local  temperature  in  the  cloud. 


the  spherical  cloud  is  too  far  away  from  the  source.  If  the  inner 
radius  is  larger  than  about  10-20  pc,  the  luminosities  of  the 
[O  iv]  and  [Ne  v]  lines  drop  dramatically.  Interestingly,  the  He  ii 
line,  though  similar  to  [O  iv]  in  ionization  potential,  is  nearly 
insensitive  to  the  inner  cloud  radius.  This  is  likely  because  the 
[O  iv]  and  [Ne  v]  are  in  a  higher  stage  of  excitation  than  He  ii  and 
as  a  result  are  much  more  affected  by  dilution  of  the  ionizing 
photons  (in  other  words  a  drop  in  the  ionization  parameter) 
from  increasing  the  inner  radius.  We  conclude  from  this  that 
increasing  the  inner  radius  will  only  worsen  our  original  results 
and  will  not  affect  the  strength  of  the  low-ionization  lines.  Thus, 
we  conclude  that  the  inner  radius  is  likely  close  to  the  source  as 
our  original  input  parameter  supposes. 

Varying  the  stellar  companion  (the  05V  star).  To  test  the 
effect  of  the  stellar  companion  on  the  IR  lines  we  ran  the  base 
model  with  an  05V  star  replacing  the  B2Ib  supergiant.  The 
temperature  and  luminosity  of  the  05V  star,  42,000  K  and  3.2  x 
1039  g-i  ^  were  chosen  to  match  the  model  used  by  KWZ.  The 
results  (the  lower-right  panel  of  Eigure  5)  show  that  the  05V 
star  has  a  very  strong  impact  on  the  high-ionization  lines.  The 
luminosity  of  the  [Oiv]  line  increases  by  a  factor  of  3,  and  the 
[Ne  v]  line  by  factor  of  2.  If  the  companion  is  in  fact  an  05  V  star, 
our  measured  optical  monitor  (OM)  and  GALEX  fluxes  from 
Paper  I  are  underluminous  and  are  likely  affected  by  reddening. 
Inclusion  of  an  05V  star  also  has  a  very  different  effect  on 
the  [O  iv]  versus  the  He  ii  line.  Indeed,  the  He  ii  luminosity  is 
predicted  to  increase  by  only  11%.  The  strong  impact  of  the  05  V 
star  can  be  explained  if  the  O  star  produces  copious  O^"^  ions 
(i.p.  35  eV),  which  are  in  turn  ionized  further  by  the  ULX.  The 
O  star  thus  acts  as  a  catalyst  to  the  higher-energy  photons  from 
the  ULX.  A  comparison  between  the  oxygen  ions  structure  for 
these  models  is  presented  in  Eigure  8(a).  The  number  of  photons 
in  the  24.6-54.4  eV  interval  is  3  times  higher  for  the  05V  star 
(7.27  X  10"^^  s“^  compared  to  2.42  x  10"^^  s“^)  than  the  B2Ib 
supergiant.  Eigure  8(b)  shows  that  the  local  cloud  temperature 
for  the  05  V  star  is  also  about  2  times  higher  in  the  region  where 
the  O^"^  ions  are  produced  (at  offsets  of  10-20  pc).  Indeed,  the 


[O  III]  A. 5007  luminosity  estimated  from  CLOUDY  increases  by 
a  factor  of  6  when  the  05V  star  is  used  instead  of  the  B2Ib 
supergiant. 

The  companion  star  therefore  has  an  indirect  but  strong 
impact  on  the  [Oiv]  line  emission.  A  similar  effect  can  be 
seen  on  a  smaller  scale  for  the  [Ne  v]  line,  which  has  an  even 
higher  ionization  potential  (97  eV).  Thus,  the  [Oiv]  line  does 
not  act  simply  as  a  high-energy  photon  counter  as  is  the  case 
with  the  Hen  line,  but  can  also  provide  information  about 
abundances  of  lower  energy  photons.  This  effect  can  be  used 
to  explain  the  detection  of  [O  iv]  in  star- forming  regions.  In  this 
scenario,  [O  iv]  can  be  produced  by  relatively  faint  sources  of 
soft  X-rays  and  UV  photons  hidden  in  the  star-forming  region, 
including  such  sources  as  ULXs,  X-ray  binaries,  or  SNR.  This 
is  a  significant  result,  since  the  origin  of  [O  iv]  in  star-forming 
galaxies  has  been  the  subject  of  debate  in  the  literature  (Lutz 
et  al.  1998;  Schaerer  &  Stasinska  1999). 

Because  our  CLOUDY  modeling  overpredicts  most  of  the 
high-ionization  line  fluxes  when  an  OV5  star  is  used,  we 
conclude  that  an  OV5  star  is  likely  not  the  companion  star  for  the 
Holmberg  II  ULX.  On  the  other  hand,  later-type  stars,  such  as  an 
08  star  (or  later),  might  better  predict  the  [O  iv]  fiux  in  the  base 
model.  In  support  of  this,  if  the  ULX  is  powered  by  a  stellar- 
mass  BH,  the  companion  required  to  reproduce  the  observed 
X-ray  luminosities  must  be  very  massive  (Podsiadlowski  et  al. 
2003).  Thus  a  stellar-mass  BH  could  explain  both  the  [Oiv] 
luminosity  and  the  very  high  X-ray  luminosity.  The  fact  that 
the  Holmberg  II  ULX  resides  just  outside  a  UV-bright  region  of 
star  formation  (see  Paper  I),  also  suggests  that  the  environment 
is  well  suited  for  formation  of  such  massive  stars  (similar  to 
NGC4559  X7,  see  Soria  2007).  However,  without  the  BH  mass 
constrained,  an  IMBH  with  a  massive  companion  is  still  also  a 
plausible  option. 

3.4.  Varying  the  Geometry 

As  mentioned  at  the  beginning  of  this  section,  KWZ  found 
evidence  from  low-ionization  optical  lines  that  the  ionized 
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Figure  9.  This  figure  illustrates  the  difference  between  a  matter-bounded  and  radiation-bounded  nebula.  The  lines  are  only  labeled  in  the  right  panel,  (a)  The  dependence 
of  the  simulated  lines  from  CLOUDY  on  the  cloud  hydrogen  column  density  for  the  base  model.  The  equivalent  depth  of  the  spherical  cloud  is  also  marked  on  the 
top  Y-axis.  Both  vertical  continuous  lines  mark  the  column  density  for  the  base  model,  but  one  value  was  obtained  by  fitting  the  X-ray  data  with  the  base  model 
(matter-bounded  nebula),  while  the  other  is  the  column  in  the  CLOUDY  simulation  where  the  local  temperature  reached  4000  K  (radiation-bounded).  We  also  marked 
with  a  vertical  dashed  line  the  limit  where  the  predictions  for  the  high- excitation  lines  are  still  compatible  with  observations,  (b)  The  simulated  lines  using  the  family 
of  broken  power-law  models  fitted  to  the  X-ray  data.  The  BPL  is  marked  by  the  vertical  line.  The  bottom  X-axis  shows  the  low-energy  photon  index  for  the  models. 
On  the  top  X-axis  we  mark  the  corresponding  log  column  densities  (lower  values)  and  cloud  thickness  (upper  values).  All  these  models  are  matter-bounded. 


nebula  is  matter-bounded,  at  least  in  the  east.  This  is  also 
suggested  by  the  IR  line  maps  presented  in  Paper  I  and  by 
the  overprediction  of  the  low-ionization  lines  by  the  radiation- 
bounded  base  SED  model  above.  Paper  I  also  shows  that  the  low- 
ionization  line  fluxes  are  likely  contaminated  by  circumstellar 
emission.  We  also  noticed  that  the  absorption  suggested  by 
the  X-ray  spectrum  is  too  low-compared  to  that  implied  by  a 
radiation-bounded  geometry.  More  precisely,  for  our  base  SED 
model,  the  simulation  requires  a  column  density  of  4.2  X  10^1 
cm“^,  which  is  more  than  10  times  larger  than  the  measured 
values  from  spectral  fits  (see  also  the  beginning  of  this  section). 
This  suggests  that  the  nebula  is  also  matter-bounded  in  the  line- 
of- sight  direction. 

To  simulate  a  matter-bounded  geometry,  we  first  ran  the 
base  SED  model  with  the  geometry  described  by  a  spherical 
cloud  of  varying  thickness,  so  that  the  total  hydrogen  column 
ranges  between  10^^  and  10^^  cm“^.  These  results  are  plotted  in 
Figure  9(a).  The  column  obtained  from  the  X-ray  fit  is  shown 
as  a  vertical  line  and  labeled  “Matter  Bounded.”  For  this  value, 
the  predicted  fluxes  for  both  [Oiv]  and  Hen  lines  are  much 
lower  than  the  observed  values.  Therefore,  this  simple  matter- 
bounded  geometry  using  the  base  SED  model  is  not  sufficient 
to  reproduce  the  observed  lines.  As  Figure  9(a)  shows,  the 
minimum  column  density  required  for  the  base  model  to  predict 
the  measured  values  for  the  high-excitation  lines  is  x  10^^, 
or  a  cloud  thickness  of  16  pc  (the  dashed  vertical  line).  This 
values  is  within  the  15-25  pc  radii  observed  in  the  He  ii  emission 
line  maps  from  KWZ.  Therefore  the  matter-bounded  case  could 
work,  if  the  observed  X-ray  column  density  were  slightly  higher 
than  the  X-ray  predicted  value  of  3.26  x  10^^  cm“^. 

We  next  replaced  the  modified  PLMCD  model  for  the  BPL 
model  as  the  X-ray  fit  for  the  SED.  The  broken  power-law 


model  is  better  suited  than  the  Modified  PLMCD  model  to 
test  the  matter-bounded  geometry.  By  varying  this  absorption 
parameter  in  the  model  within  the  errors  allowed  by  the  fit, 
we  obtained  a  family  of  broken-power  models,  all  consistent 
with  the  X-ray  data,  but  predicting  very  different  slopes  in  the 
ionizing  part  of  the  SED  (the  photon  index  varies  between  1.3 
and  2.2).  A  subset  of  the  broken  power-law  models  is  shown  in 
Figure  4.  We  re-ran  our  CLOUDY  simulation  for  these  models 
with  the  corresponding  column  densities  from  the  X-ray  fits.  The 
results  are  plotted  in  Figure  9(b).  The  best-fit  broken  power-law 
X-ray  model  (BPL)  is  shown  as  a  vertical  line  in  this  figure. 
The  predicted  line  luminosities  for  this  model  are  given  in  the 
lower-left  panel  of  Figure  5. 

As  was  the  case  with  the  original,  radiation-bounded  base 
SED  model,  the  predictions  are  very  good  for  the  high-ionization 
lines.  The  cloud  radius  is  17.7  pc,  which  is  consistent  with  the 
size  of  the  nebula  found  by  KWZ  in  the  east  and  south.  However, 
contrary  to  the  original  radiation-bounded  base  model,  the  ULX 
makes  no  contribution  to  the  low-ionization  lines.  This  is  true 
even  for  the  largest  column  density,  9.2  x  10^^  cm  ^  (a  29.7  pc 
cloud  radius  for  a  gas  density  of  10  cm“^).  However,  the  spectral 
maps  in  Figure  3  of  Paper  I  show  that  the  ULX  is  likely  to 
contribute  at  least  to  [Ne  iii]  and  probably  [Si  ii] . 

Thus,  we  conclude  that  the  matter-bounded  geometry  alone 
does  not  fit  all  of  the  mid-IR  spectral  observations  well.  In  fact, 
these  results  suggest  that  the  geometry  is  complicated,  being 
asymmetric  and  including  both  radiation-  and  matter-bounded 
geometries.  KWZ  found  similar  results.  They  find  that  to  the 
south  and  east.  He  ii  and  H^  emission  suggests  a  matter-bounded 
geometry:  the  low-excitation  lines  are  absent  here  but  detected 
in  the  west  at  larger  offsets.  Our  CUBISM  maps  from  Paper  I 
suggest  something  similar,  though  the  spatial  resolution  here 
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Figure  10.  Mappings  III  shocks  and  photoionization  models.  The  shock  grids  are  plotted  for  magnetic  parameters:  0.5,  1,  2,  3.23,  5,  and  10  /xG  cm^/^  and  for 
shock  velocities  125,  250,  300,  400,  500,  700,  and  1000  km  s“^  (a)  [Oiv]  line  intensity  shocks  predictions.  The  shocks  plus  precursor  models  vary  little  with  the 
magnetic  parameter,  and  we  only  show  the  model  with  a  nominal  equipartition  value  (3.23  /xG  cm^/^,  the  dashed  line).  The  horizontal  lines  correspond  to  the  brightest 
(continuous)  and  respectively  faintest  (dotted)  pixels  in  the  [O  iv]  line  map  in  Figure  3  of  Paper  I.  (b)  Optical  lines  diagnostic  with  shocks  and  photoionization  models. 
The  pure  shock  models  are  in  red,  shocks  plus  precursor  in  blue,  and  the  photoionization  grid  in  black.  For  the  photoionization  grid  we  varied  the  density  (1,  10,  and 
100  cm“^),  and  the  inner  disk  temperature  (0.1,  0.2,  0.4,  0.8,  and  1.0  keV),  and  assumed  a  metallicity  of  0.1  Zq.  For  Holmberg  IX  ULX  and  MF16  we  show  with 
arrows  the  approximate  locations  for  Solar  metallicity.  While  Holmberg  IX  ULX  is  consistent  with  shocks  with  modest  velocities  (<100  km  s“^),  Holmberg  II  ULX 
and  MF16  seem  to  require  very  high  velocities.  They  are  however  consistent  with  photoionization. 


is  limited  compared  with  the  HST  observation.  The  fact  that 
the  low-ionization  line  luminosities  predicted  by  the  matter- 
bounded  models  are  much  lower  than  those  predicted  by  the 
radiation-bounded  geometry  suggests  that  a  combination  of  the 
two  geometries  would  yield  predicted  values  that  are  very  close 
to  the  observations. 

4.  SHOCK  MODELS 

The  radiation-bounded/matter-bounded  geometric  scenario 
presented  in  Section  3  to  explain  the  mid-IR  line  fluxes  in  the 
Holmberg  II  ULX  presumes  photoionization  as  the  source  of  the 
line  emission.  However,  shocks  could  also  produce  such  lines. 
Indeed,  a  combination  of  X-ray  photoionization  and  shocks  was 
recently  proposed  to  explain  the  morphology  of  optical  lines 
observed  in  the  vicinity  of  EMC  X-1  (Cooke  et  al.  2007). 
In  this  case  the  shocks  are  likely  driven  by  a  jet.  Abolmasov 
et  al.  (2007)  performed  optical  spectroscopy  on  eight  ULXs 
with  nebular  counterparts,  including  the  Holmberg  II  ULX, 
and  found  evidence  for  shocks  in  all  cases  based  on  the  [S  ii] 
A.A.67 17,673 1/Ha  ratio.  However,  using  high-ionization  lines 
such  as  He  ii  and  [O  iii],  these  authors  also  found  that  Holmberg 
II  ULX  and  three  other  ULXs  required  photoionization  to 
sufficiently  explain  the  optical  ratios. 

Lehmann  et  al.  (2005)  found  evidence  of  radial  velocity 
variations  of  ±50  km  s~^  from  a  region  similar  of  size  to  the 
He  III  nebula.  It  is  therefore  possible  that  part  of  the  ionization 
is  caused  by  shocks.  Lutz  et  al.  (1998)  suggested  that  velocity 
shocks  of  100-200  km  s“^  can  produce  significant  [Oiv] 
emission.  However,  when  we  also  take  into  account  the  optical 
data,  and  in  particular  the  detection  of  the  He  ii  recombination 
line,  our  data  provide  strong  evidence  against  shocks  as  the 
source  of  the  [Oiv]  emission.  This  is  primarily  because  the 
He  II  is  not  easily  produced  by  shocks.  Instead,  the  simultaneous 
detection  of  these  lines  argues  for  the  presence  of  significant 


radiation  over  the  Hen  edge  (see  also  Schaerer  &  Stasihska 
1999),  and  therefore  suggests  photoionization  as  the  most  likely 
line  production  mechanism. 

Nevertheless,  we  checked  to  see  if  shock  models  could  explain 
the  mid-IR  line  luminosities  and  ratios  observed  in  Holmberg  II. 
We  used  Mappings  III  grids  from  Allen  et  al.  (2008)  to  construct 
diagnostic  diagrams.  The  predicted  [O  iv]  emission  from  shocks 
is  shown  in  Ligure  10(a).  We  show  both  pure  shocks  and 
shocks  with  ionized  precursors  at  a  gas  density  of  1  cm“^, 
and  metallicity  0.1  Zq.  Lor  these  models,  the  parameters  varied 
are  the  shock  velocity  and  the  “magnetic  parameter” 
where  B  is  the  magnetic  field  and  n  is  the  gas  density  (Dopita  & 
Sutherland  1996). 

The  CUBISM  measured  fluxes  from  the  brightest  and  faintest 
pixels  in  the  [Oiv]  map  (yellow  and  blue  pixels  in  Ligure  3 
of  Paper  I,  respectively)  are  shown  as  horizontal  lines  in 
Ligure  10(a).  We  note  that  the  maximum  intensity  of  the  line 
is  likely  higher  than  these  levels.  Indeed  the  [O  iv]  emission  is 
not  spatially  resolved  and  the  plotted  levels  are  actually  averages 
over  the  size  of  the  LH  pixel  (67.5  pc),  rather  than  the  maximum 
from  each  of  the  averaged  slits.  In  spite  of  this,  the  levels  plotted 
in  Ligure  10(a)  requires  significant  shock  velocities  and  such 
velocities  (>200  km  s“^)  were  not  detected  by  Lehmann  et  al. 
(2005). 

Lor  optical  lines  we  chose  the  He  ii  and  the  [O  iii] 
Ak4363,5007  lines,  as  per  Evans  et  al.  (1999).  The  ratio  [Oiii] 
A.4363/[Oiii]  A.5007  is  sensitive  to  the  gas  temperature.  Using 
the  data  from  Lehmann  et  al.  (2005),  we  calculated  the  ratio 
([Oiii]  14959  [Oiii]  15007)/[Oiii]  14363  =  50,  and  esti¬ 
mate  a  temperature  of  17,800  K  (Osterbrock  &  Lerland  2006). 
This  is  consistent  with  the  temperature  of  the  gas  in  the  region 
of  the  cloud  where  O^"^  ions  are  produced  (see  Ligure  8(b)). 
Lor  the  Heii/H^  ratio  we  chose  an  average  (0.2)  between  the 
measurements  of  KWZ,  Pakull  &  Mirioni  (2002),  and  Lehmann 
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et  al.  (2005).  The  shock  grids  for  the  optical  lines  are  plotted 
in  Figure  10(b).  For  comparison,  we  also  plot  a  photoionization 
grid  obtained  with  Mappings  III,  using  parameters  similar  to 
those  used  in  the  CFOUDY  modeling.  Finally,  we  plot  previ¬ 
ously  published  data  for  two  other  UFXs  with  detected  nebulae, 
Holmberg  IX  UFX  and  MF16  in  NGC  6946  (Abolmasov  et  al. 
2007).  While  the  Holmberg  IX  UFX  data  are  consistent  with 
shocks  with  modest  velocities  (<  100  km  s- ' ),  Holmberg  II ULX 
and  MF16  seem  to  require  very  high  velocities.  Even  for  the 
shock  plus  precursor  models,  velocities  in  excess  of  300  km  s“^ 
are  required  to  reproduce  the  observed  line  ratios  and  Fehmann 
et  al.  (2005)  found  velocities  of  only  50  km  s“^ 

In  summary,  some  contribution  from  shocks  to  the  high- 
excitation  lines  detected  from  Holmberg  II  UFX  cannot  be  ruled 
out  completely.  However,  if  there  is  some  contribution  from 
shocks  to  the  high-ionization  emission  line  luminosities,  it  is 
likely  very  small,  since  both  He  ii  and  [O  iv]  lines  require  shocks 
with  velocities  in  excess  of  200  km  s“^ ,  and  such  velocities  were 
not  seen  by  Fehmann  et  al.  (2005). 

5.  DISCUSSION  OF  THE  MODEFING  RESUFTS 

We  find  that  the  Spitzer  observations  of  the  Holmberg  II 
UFX  and  especially  the  detection  of  the  [O  iv]  line  is  consistent 
with  photoionization  by  radiation  from  the  UFX.  Of  all  the 
lines  examined,  only  the  [Oiv]  was  shown  in  Paper  I  to  be 
correlated  with  the  UFX.  However,  we  also  found  in  Paper  I 
that  some  contribution  by  the  UFX  to  the  [Sin],  [Siii],  and 
[Neiii]  emission  line  flux  is  also  likely.  The  modeling  results 
presented  here  are  based  on  the  observations  presented  in  Paper  I 
and  are  summarized  as  follows. 

1.  The  base  model  (modified  PLMCD).  The  CFOUDY  sim¬ 
ulations  show  that  the  Spitzer  spectrum  is  consistent  with 
photoionization  by  the  accretion  powered  emission  from  the 
UFX.  We  find  that  limited  sensitivity  might  explain  the  rel¬ 
atively  low  fiuxes  observed  for  some  low-ionization  lines. 
The  base  SED  model  predicts  a  slightly  lower  [O  iv]  lumi¬ 
nosity  than  the  measured  values  for  the  radiation-bounded 
case  (Figure  5),  and  even  lower  for  the  matter-bounded 
case  (Figure  9(a)).  Previous  X-ray  results  have  shown  that 
the  accretion  disk  temperature  is  likely  lower  than  in  our 
base  model  (KWZ  obtained  "^0.2  keV).  Such  a  disk  temper¬ 
ature  will  increase  the  predicted  line  luminosity  to  match 
our  measurements  (Figure  6(a)).  A  slightly  higher  metal- 
licity  (0.2  Zq)  can  yield  better  fits  to  both  the  [Oiv]  and 
He  II  data  for  this  base  model.  This  is  because  the  He  ii  line 
is  already  overpredicted  and  raising  the  metallicity  will  not 
increase  this  discrepancy  (Figure  7(a)).  In  support  of  this,  a 
metallicity  of  0.17  was  measured  by  Pilyugin  et  al.  (2004). 
We  note  that  while  this  model  reproduces  the  spectra  very 
well,  it  requires  absorbing  columns  in  the  X-rays  that  are 
inconsistent  with  the  observations. 

2.  The  PLMCD  model.  We  find  that  the  widely  used  PFMCD 
model,  when  extrapolated  to  the  UV,  is  not  consistent  with 
the  observations.  The  high-ionization  lines  predicted  using 
this  model  are  at  least  10  times  larger  than  observed.  In 
this  case,  lowering  the  metallicity  to  match  the  [O  iv]  will 
not  work,  because  Hen  will  still  be  much  overpredicted. 
This  model  has  been  used  extensively  to  fit  the  X-ray 
data  from  UFXs  (e.g.,  “cool  disks,”  Miller  et  al.  2004). 
However,  the  X-ray  community  has  more  recently  begun 
to  fit  these  spectra  with  more  physically  motivated  models 
(e.g.,  Freeland  et  al.  2006;  Done  &  Kubota  2006;  Kuncic 


et  al.  2007 ;  Soria  2007 ;  Soria  &  Kuncic  2008),  which  appear 
to  be  consistent  with  high-quality  UFX  X-ray  data  (Berghea 
et  al.  2008).  That  the  PFMCD  model  is  inconsistent  with 
the  Spitzer  observations  of  the  Holmberg  II  UFX  supports 
these  previous  X-ray  studies. 

3.  Matter  versus  radiation-bounded  geometry.  The  intensity 
and  morphology  of  the  IR  lines  provides  spatial  information 
about  the  emission  line  nebula  around  the  Holmberg  II 
UFX.  The  geometry  consistent  with  our  infrared  data 
(and  previously  published  optical  data)  is  approximately 
spherical  but  not  symmetrical.  The  CUBISM  maps  show 
that  to  the  east  and  south  the  nebula  shows  only  high- 
excitation  lines  that  are  likely  matter-bounded,  while  in 
the  west,  the  detection  of  lower  excitation  lines  might 
suggest  a  radiation-bounded  cloud.  However,  observational 
effects,  such  as  aperture  and  sensitivity  effects,  preclude 
any  definitive  conclusion  on  this  model.  The  discrepancies 
found  in  Section  3.1  between  the  CFOUDY  model  and  the 
measured  values  for  the  low-excitation  lines  suggest  that 
the  matter-bounded  geometry  is  at  play  at  least  in  part 
(probably  to  east  and  south  as  found  from  the  spectral 
maps).  Moreover,  the  absorption  measurements  from  the 
X-ray  data  require  the  cloud  to  be  at  least  matter-bounded 
toward  our  line  of  sight.  Therefore,  we  find  that  both  matter- 
bounded  and  radiation-bounded  geometries  are  likely  at 
work  to  describe  the  observed  spectra  in  both  the  Mid-IR 
and  X-rays. 

4.  Shock  models.  The  shocks  models  for  both  IR  and  optical 
emission  lines  show  that  shocks  with  velocity  >200  km  s“^ 
are  required  to  generate  both  the  Hen  and  [Oiv]  lines 
observed  in  Holmberg  II.  The  detected  radial  velocity 
variations  around  the  UFX  are  of  the  order  of  ±50  km  s“^ 
(Fehmann  et  al.  2005).  Contribution  from  shocks  to  these 
lines  is  therefore  likely  to  be  small. 

5.  True  X-ray  luminosity.  The  UFX  luminosity  estimated  us¬ 
ing  photoionization  modeling  is  independent  of  the  estimate 
based  on  the  X-ray  flux,  which  can  be  affected  by  both  ab¬ 
sorption  and  beaming  (King  et  al.  2001).  The  nebular  emis¬ 
sion  line  luminosities,  on  the  other  hand,  are  time-averaged 
true  luminosities  and  can  provide  information  about  the 
geometry  of  emission  and  the  surrounding  ISM. 

It  is  now  widely  accepted  that  UFXs  are  variable  on 
timescales  from  days  to  years.  These  timescales  are  much 
shorter  than  the  ionization  equilibrium  timescales  in  the 
nebula.  Chiang  &  Rappaport  (1996)  investigated  time- 
dependent  photoionization  of  H  and  He  for  supersoft  X-ray 
sources.  The  Hen  ionization  equilibrium  at  the  90%  level 
is  reached  in  ^^0.7  Tr  after  source  tum-on,  where  Tr  = 
3000  yr  is  the  recombination  time  for  He^"^  from  KWZ.  For 
periodic  sources,  it  was  found  that  the  He  ii  line  luminosity 
decreases  significantly  with  increasing  source  period.  For 
example,  for  a  period  of  10  and  a  duty  cycle  of  10%, 
the  estimated  luminosity  is  4.7%  of  that  expected  from  a 
steady  source  with  the  same  peak  luminosity.  That  is  a 
53%  deviation  from  a  linear  scaling  law  (10%  of  the  steady 
state).  This  implies  that  if  the  Holmberg  II  UFX  is  variable, 
the  estimates  of  KWZ  based  on  the  He  ii  line  are  actually 
lower  limits  of  the  time-averaged  luminosity. 

For  the  recombination  time  is  shorter,  Tr  ~  100  yr 
(Osterbrock  &  Ferland  2006),  for  the  same  ionization 
conditions:  electron  density  10  cm“^  and  temperature 
20,000  K.  We  expect  therefore  that  the  [O  iv]  line  reaches 
equilibrium  in  a  shorter  time  (<100  yr),  but  still  long 
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enough  to  provide  a  time- averaged  estimate  of  the  true 
luminosity  of  the  ULX  over  ~10-100  yr.  We  note  that 
light  travel  time  (53  yr  for  a  nebula  of  16  pc)  is  an 
important  factor  in  this  case,  because  it  is  comparable  to 
the  recombination  time. 

Einally,  the  morphology  of  the  emission  lines,  from  both 
the  previously  published  optical  data  and  our  IR  data  sug¬ 
gest  emission  that  is  nearly  isotropic,  being  inconsistent 
with  narrow  beaming.  Moreover,  photoionization  model¬ 
ing  with  our  base  model  is  consistent  with  the  detected 
high-ionization  lines,  the  true  bolometric  luminosity  of  this 
model  is  >10^^  erg  s“^  This  is  well  above  the  Eddington 
luminosity  for  a  stellar-mass  BH.  In  addition,  some  obser¬ 
vational  results  indicate  that  the  BH  in  Holmberg  II  is  likely 
much  larger  than  a  stellar-mass  BH  (see  the  conclusion  be¬ 
low). 

6.  Impact  on  IR  studies  ofstarbursts.  We  find  that  the  predicted 
[O  iv]  line  luminosity  is  significantly  affected  by  the  type 
of  the  companion  star.  Eaint  [Oiv]  has  been  detected  in 
many  starbursts  and  star-forming  regions,  but  its  origin  is 
still  under  debate.  Our  analysis  provides  a  new  mechanism 
for  this  emission,  suggesting  that  the  [O  iv]  line  could  be 
produced  by  relatively  faint  X-ray  sources  including  X-ray 
binaries,  SNR,  or  ULXs  in  starburst  galaxies. 

6.  CONCLUSION 

The  bolometric  luminosity  of  the  base  model  is  1.33  x 
10"^^  erg  s“^.  Eor  accretion  within  the  Eddington  limit,  the 
mass  of  the  central  BH  is  at  least  85  Mq.  The  estimate  based 
on  the  inner  disk  temperature  (0.38  keV)  in  the  MCD  model 
gives  a  much  larger  mass,  994  Mq.  Einally,  scaling  from  AGN 
estimates,  the  [Oiv]  emission  predicts  a  mass  as  large  as 
10"^  Mq.  All  of  the  estimates  presented  here  presume  the  BH 
is  accreting  within  the  Eddington  limit.  However,  while  our 
analysis  excludes  strong  beaming  in  Holmberg  II  ULX,  a  sub- 
Eddington  IMBH  or  a  super-Eddington  stellar-mass  BH  are  both 
plausible  ionization  mechanisms  that  stem  from  our  analysis. 
Indeed  we  show  in  Section  3.3  that  a  massive  companion  (e.g., 
B2Ib  star)  and  a  stellar-mass  BH  can  reproduce  the  [Oiv] 
luminosities  see  in  the  Spitzer  spectra,  as  can  an  IMBH  and 
a  similar  companion.  The  CLOUDY  modeling  results  alone 
cannot  discriminate  between  these  two  sources  or  the  BH  mass 
estimates,  since  we  have  no  way  of  constraining  the  accretion 
rate. 
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